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Abstract 

The determination of the ordering of the neutrino masses (the hierarchy) is probably a crucial 
prerequisite to understand the origin of lepton masses and mixings and to establish their relation- 
ship to the analogous properties in the quark sector. Here, we follow an alternative strategy to 
the usual neutrino-antineutrino comparison in long baseline neutrino oscillation experiments: we 
exploit the combination of the neutrino-only data from the NOz^A and the T2K experiments by 
performing these two off-axis experiments at different distances but at the same {E)/L, where {E) 
is the mean neutrino energy and L is the baseline. This would require a minor adjustment to the 
proposed off-axis angle for one or both of the proposed experiments. 

PACS numbers: 14.60Pq 
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I. INTRODUCTION 



During the last several years the physics of neutrinos has achieved a remarkable progress. 
The experiments with solar [1-6], atmospheric [7], reactor [8], and also long-basehne ac- 
celerator [9] neutrinos (improved with the recent MINOS [10] results [11]) have provided 
compelling evidence for the existence of neutrino oscillations, implying non zero neutrino 
masses. The present data^ require two large [612 and ^23) and one small (^^13) angles in the 
neutrino mixing matrix [14], and at least two mass squared differences, Am^j = rrtj — mf 
(where m^'s are the neutrino masses), one driving the atmospheric {Am^i) and the other 
one the solar (Am|i) neutrino oscillations. The mixing angles Ou and 623 control the solar 
and the dominant atmospheric neutrino oscillations, while ^13 is the angle limited by the 
data from the CHOOZ and Palo Verde reactor experiments [15, 16]. 

The Super-Kamiokande (SK) [7] and K2K [9] data arc well described in terms of dominant 
~^ {i^n ~^ ^t) vacuum oscillations. The MINOS Collaboration has recently reported the 
first neutrino oscillation results from 1.27 x 10^° protons on target exposure of the MINOS 
far detector [11]. The value of the oscillation parameters from MINOS are consistent with 
the ones from K2K, as well as from SK data. A recent global fit [17] (see also Ref. [18]) 
provides the following 3(7 allowed ranges for the atmospheric mixing parameters 

lAmiil = (1.9 - 3.2) X 10-^eV^ 0.34 < sin^ ^23 < 0.68 . (1) 

The sign of Amg^, sign(Am3i), cannot be determined with the existing data. The two 
possibilities, Am|]^ > or Amg^^ < 0, correspond to two different types of neutrino mass 
ordering: normal hierarchy and inverted hierarchy. Both possibilities are illustrated in 
Fig. 1, extracted from Ref. [19]. In addition, information on the octant in which ^23 lies, if 
sin^2^^23 7^ 1, is beyond the reach of present experiments. 

The 2-neutrino oscillation analysis of the solar neutrino data, including the results from 
the complete salt phase of the Sudbury Neutrino Observatory (SNO) experiment [6], in 
combination with the KamLAND spectrum data [20], shows that the solar neutrino oscilla- 

^ We restrict ourselves to a three-family neutrino analysis. The unconfirmed LSND signal [12] cannot be 
explained within this context and might require additional light sterile neutrinos or more exotic explana- 
tions. The ongoing MiniBooNE experiment [13] is going to test the oscillation explanation of the LSND 
result. 
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tion parameters lie in the low-LMA (Large Mixing Angle) region, with best fit values [17] 
Amii = 7.9 X 10-^ eV^ and sin^ 9^2 = 0.30. 

A combined 3-neutrino oscillation analysis of the solar, atmospheric, reactor and long- 
bascline neutrino data [17] constrains the third mixing angle to be sin^ ^^13 < 0.041 at the 
3(7 C.L. However, the bound on sin^ 6*13 is dependent on the precise value of Amg^. 

The future goals for the study of neutrino properties is to precisely determine the already 
measured oscillation parameters and to obtain information on the unknown ones: namely 
^13, the CP-violating phase 5 and the type of neutrino mass hierarchy (or equivalently 
sign(Amii)). 

In this paper we concentrate on the extraction of the neutrino mass hierarchy by combin- 
ing the Phase I (neutrino-data only) of the long-baselinc appearance experiments T2K [21] 
and NOi/A [22], both exploiting the off-axis technique^. We start by presenting the general 
formalism in Sec. II. In Sec. Ill we describe the strategy and the experimental setups. In 
Sec. IV the experimental configuration is optimized in order to maximize the sensitivity to 
the mass hierarchy extraction and we present our results for different possible experimental 
combinations. Finally, in Sec. V, we make our final remarks. For our analysis, unless other- 
wise stated, we will use a representative value of lAmg^l = 2.4 x 10^^ eV^ and sin^ 2^^23 = 1- 
For the solar oscillation parameters Am|i and ^12, we will use the best fit values quoted 
earlier in this section. 

II. THE OSCILLATION PROBABILITY 

The mixing angle ^13 controls v^j, v^, and v^j, v^, conversions in long-baseline v^, 
appearance experiments and the v^. disappearance in short-baseline reactor experiments. 
Present and future conventional beams [10], super-beams with an upgraded proton source 
and a more intense neutrino flux and future reactor neutrino experiments have the possibility 

^ A neutrino beam with narrow energy spectrum can be produced by placing the detector off-axis, i. e., 
at some angle with respect to the forward direction. The resulting neutrino spectrum is very narrow in 
energy (nearly monochromatic, A.E/E ^ 15 — 25%) and peaked at lower energies with respect to the 
on-axis one. The off-axis technique allows a discrimination between the peaked oscillation signal and 
the intrinsic background which has a broad energy spectrum [23] . 
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Fractional Flavor Content varying cos 6 



FIG. 1: The range of probability of finding the a-flavor in the i-th mass eigenstate as indicated 
for the two different mass hierarchies for the best fit values of the solar and atmospheric mixing 
parameters, sin^ 2^13 = 0.1 as the CP-violating phase, 6 is varied. 

to measure^, or set a stronger limit on, ^13. Therefore, with the possibility of the first 
measurement of ^13 being made by a 1-to 2-km baseline reactor experiment [30], the long- 
baseline off-axis z/e appearance experiments, T2K [21] and NOz/A [22], need to adjust their 
focus to emphasize other physics topics. The most important of these questions is the form 
of the mass hierarchy, normal (Amg^^ > 0) versus inverted (Am|^ < 0), and whether or not 
leptonic CP violation occurs. 

Consider the probability -P(z^^ — > i^e) in the context of three-neutrino mixing in the 
presence of matter [31], represented by the matter parameter a, defined as a = Gf^b/v^, 
where is the average electron number density over the baseline, taken to be constant 
throughout the present study. Defining Ajj = ^ convenient and precise approximation 

is obtained by expanding to second order in the following small parameters: ^13, A21/A32, 
A21/ aL and A21. The result is (details of the calculation can be found in Ref. [25], see also 

^ Smaller values of this mixing angle could be accessed by very long baseline experiments such as neutrino 
factories [24, 25], or by /3-beams [26-29] which exploit neutrinos from boosted-ion decays. 
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Ref. [32]): 



sin ^23 sin 2^13 ( ^l ) ^^^^^^i - aL)e-^(^32+'^) + cos^23 sin 2^12 ] sin (aL) 



= sin^ 



^23 sin^ 2^13 (^^^T^) sin' (A31 - aL) 



+ sin 2^13 sin 2^23 sin 2^12 f — ) ( . ^'^^ ^ ] sin (aL) sin (A31 — aL) cos (A32 + 5) 

\aL J VA31 - aL ' 



+ cos" ^23 Sin' 2^12 y-^j sin' (aL) , (2) 

where L is the basehne and a —a, 6 —S for P{i^ij, —>■ t'e). 

Suppose — > i^e) < ^ J^e) for given energy and basehne, then in vacuum this 

imphes CP violation. However, in matter, this imphes CP violation only for the normal 
hierarchy but not necessarily for the inverted hierarchy around the first oscillation maxi- 
mum. The different index of refraction for neutrinos and antineutrinos induces differences 
in the u, v propagation that could be misinterpreted as CP violation [33]. Typically, the 
proposed long baseline neutrino oscillation experiments have a single far detector and plan 
to run with the beam in two different modes, neutrinos and antineutrinos. In principle, 
by comparing the probability of neutrino and antineutrino flavor conversion, the values of 
the CP-violating phase b and of sign(Am|i) could be extracted. However, different sets 
of values of CP-conserving and violating parameters, (^^13, ^^23, b-, sign(Am3^)), lead to the 
same probabilities of neutrino and antineutrino conversion and provide a good description 
of the data at the same confidence level. This problem is known as the problem of degen- 
eracies in the neutrino parameter space [34-38] and severely affects the sensitivities to these 
parameters in future long-baseline experiments. Many strategies have been advocated to 
resolve this issue. Some of the degeneracies might be eliminated with sufficient energy or 
baseline spectral information [25] . In practice, statistical errors and realistic efficiencies and 
backgrounds limit considerably the capabilities of this method. Another detector [35, 39-42] 
or the combination with another experiment [43-51] would, thus, be necessary"'. 

The use of only a neutrino beam could help in resolving the type of hierarchy when 



^ New approaches which exploit other neutrino osciUations channels such as muon neutrino disappearance 
have been proposed [52] for determining the type of hierarchy. They require very precise measurements 
of the neutrino oscillation disappearance oscillation parameters though. 
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two different long-baselines are considered [42-44, 53, 54]. It was shown in ref. [44] that 
if the (E) I L for the two different experiments is approximately the same then the allowed 
regions for the two hierarchies are disconnected and thus this method for determining the 
hierarchy is free of degeneracies.^ Naively, we can understand this method in the following 
way for sin^ 2^13 > 0.01: assume that matter effects are negligible for the short baseline, 
then at the same {E)/L, if the oscillation probability at the long baseline is larger than the 
oscillation probability at the short baseline, one can conclude that the hierarchy is normal, 
since matter effects enhance the neutrino oscillation probabilities for the normal hierarchy. 
For the inverted hierarchy the oscillation probability for the long baseline is suppressed 
relative to the short baseline. This will be explained in some detail in the next section. 

III. OUR STRATEGY: ONLY NEUTRINO RUNNING AND TWO DETECTORS 

Following the hue of thought developed by Minakata, Nunokawa and Parke [44], we 
exploit the neutrino data from two experiments at different distances and at different off- 
axis locations. The off-axis location of the detectors and the baseline must be chosen such 
that the (E) /Lis the same for the two experiments. Here we explain the advantages of such 
a strategy versus the commonly exploited ncutrino-antineutrino comparison. 

Suppose we compute the oscillation probabilities Piv^ — > Ue) and P{p^ — > u^) for a given 
set of oscillation parameters and the CP-phase 5 is varied between and 27r: we obtain a 
closed CP trajectory (an ellipse) in the bi-probability space of neutrino and antineutrino 
conversion [36] . Matter effects are responsible for the departure of the center of the eUipses 
from the diagonal line in the bi-probability plane for the both cases of normal and inverted 
hierarchy, sec Fig. 2 (a), where we have illustrated the case for E = 2.0 GcV and L = 810 
km. The distance between the center of the ellipse for the normal hierarchy (lower blue) 
and that for the inverted hierarchy (upper red) is governed by the size of the matter effects. 
Notice that the eUipses overlap for a significant fraction of values of the CP-phase 5 for every 
allowed value of sin^ 2^13. This indicates that, generically, a measurement of the probabihty 
of conversion for neutrinos and antineutrinos cannot uniquely determine the type of hierarchy 
in a single experiment. This makes the determination of sign(Am3j^) extremely difficult, i. e., 

^ This method, while optimal for extracting the neutrino mass hierarchy, is not the most powerful one for 
extracting the CP violating phase 6 for which antineutrino running would be, in general, necessary. 
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FIG. 2: (a) Bi-prohahility neutrino -antineutrino ellipses at the far distances (810 km) for normal 
(lower blue) and inverted (upper red) hierarchies, (h) Same as (a) but in the bi-neutrino plane at 
the short (295 km) and far distances (810 km). 

the sign(Am3^)-extraction is not free of degeneracies. 

In the case of bi-probability plots for neutrino-neutrino modes at different distances 
(which will be referred as near (N) and far (F)), the CP-trajectory is also elliptical (see 
Fig. 2 (b)). The overlap of the two bands, which implies the presence of a degeneracy of the 
type of hierarchy with other parameters, is controlled by the difference in the slopes and the 
width of the bands. Using the fact that matter effects are small {aL ^ A31), we can perform 
a perturbative expansion and assuming that the {E) / L of the near and far experiments is 
the same^, at first order, the ratio of the slopes reads [44] 

— ~ 1 + 4(aNi^N - OfIvf) ( ; — lir-^] , (3) 

a- V^3i tan(A3i)y 

where a+ and «_ are the slopes of the center of the ellipses as one varies 6'i3 for normal and 

inverted hierarchies, and a-p and on are the matter parameters, L-p and are the baselines 

for the two experiments^. The separation between the center of the ellipses for the two 



^ The reason for this choice of {E) / L is explained in the next paragraph. 

^ Notice that although we are using the constant density approximation, ap and are different because 
the average density depends on the baseline. For T2K (NOi^A ) we use an average density times electron 
fraction equal to 1.15 (1-40) g • cm~'^. 
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hierarchies increases as the difference in the matter parameter times the path length for the 
two experiments increases. Also, since ( — cot A) is a monotonically increasing function, 
we conclude that the smaller the energy, the larger the ratio of slopes, at least for the same 

{E)IL. 

However the width of the ellipses is crucial: even when the separation between the central 
axes of the two regions is substantial if the ellipses for the normal and inverted hierarchy 
overlap the hierarchy cannot be resolved for values of CP phase, 5, for which there is overlap. 
The width of the ellipses is controlled by the difference in the {E) /L of the two experiments. 
For fixed ^13 the ellipses are fiat if there exists two different 5 which give the same oscillation 
probability for both the near, P/v, and far, Pp, detectors. That is, if 

Piv(^i3, 5) = PAr(^i3, 5') and (^13, 5) = Pf(^i3, 5') (4) 

has non-trivial solutions for 5 and 5'. Eqs. (4) are satisfied when cos(A -\- 5) — cos(A + 5') 
for both the near and far experiments simultaneously. Both equations can only be satisfied 
when Aiv = A^, that is the two experiments have the same {E)/L. 

Thus, we have two conditions to satisfy to optimize the determination of the neutrino 
mass hierarchy: 

• (a) maximize the difference in the factor (oL) for both experiments and 

• (b) minimize the ellipses width by performing the two experiments at the same {E) / L. 

The most promising way to optimize the sensitivity to the hierarchy with relatively near 
term data is therefore to focus on the neutrino running mode and to exploit the Phase I 
data of the long-baseline off-axis i/g appearance experiments, T2K and NOi/A. T2K utilizes a 
steerable neutrino beam from J- PARC and Super-Kamiokande and maybe eventually Hyper- 
Kamiokande as the far detector. The beam will peak at 0.65 GeV with the detector off-axis 
by an angle of 2.5° at 295 km. For this configuration the matter effects are small but 
not negligible [55]. NOz/A proposes to use the Fermilab NuMI beam with a baseline of 
810 km with a 30 kton low density tracking calorimeter with an efficiency of 24%. Such 
a detector would be located 12 km off-axis distance from the beam center ai L — 810 km 
(corresponding to 0.85° off-axis angle), resulting in a mean neutrino energy of 2.0 GeV. 
Matter effects are quite significant for NOi/A. Therefore, the condition (a) is satisfied, since 
(aL)NOi/A — 3(aL)T2K ■ What about the condition (b)? A back-of-the-envelope calculation 
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indicates us that the current off-axis angles are not such that (E) /L of the two experiments 
will be the same. However, by placing the detector (s) at slightly different off-axis angle(s), 
one can arrange that the (E) /L of the two experiments to be exactly the same. This strategy 
would only need half of the time of data taking (because we avoid the antineutrino running), 
when compared to the standard one (i.e. running in neutrinos and antineutrinos at a fixed 
energy, (E), and baseline, L). 

IV. OPTIMIZING THE NOi/A AND T2K DETECTOR LOCATIONS 

In this section we present what could be achieved if NOi/A and T2K setups are carefully 
chosen, focusing on the physics potential of the combination of their future data. We define 
the Phase I of the experiments as follows. For the T2K experiment, we consider 5 years 
of neutrino running and SK as the far detector with a fiducial mass of 22.5 kton and 70% 
detection efficiencies. For the NOi/A experiment, we assume 6.5 x 10^° protons on target 
per year, 5 years of neutrino running and the detector described in the previous section. 

We present in Figs. 3 the bi-event neutrino-neutrino plots, fixing the off-axis location 
of the T2K detector to 2.5° but varying the off-axis distance of the NOi^A far detector, 
where the number of events is computed by integrating the neutrino fiux convoluted with 
the cross section, the oscillation probability and the efficiency over an energy window of 1 
GeV centered around the mean neutrino energy (for both NOi^A and T2K experiments). If 
the location of the NOi^A far detector is the proposed one, i.e. 12 km off-axis, the width 
of the ellipses is not negligible and the determination of the hierarchy would become a 
really difficult task even if nature has been kind in its choice of 5, see Fig. 3 (a). As we 
anticipate in the previous section, to maximize the difference in the product aL, see Eq. (3), 
does not ensure a degeneracy- free hierarchy measurement, due to the impact of the ellipses 
width. Notice, however, that the width of the ellipses is minimal when the (E) /L of both 
experiments is almost the same: this occurs when the NOi/A far detector is located 14 
km off axis, see Fig. 3 (c), for which the mean energy (E) ~ 1.7 GeV. If the NOi/A far 
detector is placed 16 km off-axis, the {E)/L of the T2K and NOz/A experiments is no longer 
the same and the ellipses width has grown with respect to the 14 km off-axis setup. The 
most important lesson we have learned is that if the T2K beam is off- axis by an angle of 
2.5°; the ideal location for the NOvA detector to optimize the hierarchy extraction would he 
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u: 810 km 14 km (1,7 GeV) v. 810 km 16 km (1.4 GeV) 



FIG. 3: (a)Bi-event neutrino-neutrino ellipses at the NOuA and T2K experiments for normal 
(lower blue) and inverted (upper red) hierarchies. From bottom up, the dashed ellipses correspond 
to sin^ 2013 varying from 0.01 to 0.2 with a stepsize o/O.Ol and the solid ellipses illustrate the case 
sin^ 2^13 = 0.1. The T2K far detector is located off-axis by an angle of 2.5° . The NOvA detector 
is placed 12 km off-axis. (b,c,d) Same as (a) but with the NOvA detector located at 13 km, 14 km 
and 16 km off- axis, respectively. 

14 km off-axis. Fig. 3 shows the potential of the combination of the data from the Phase 
I (only neutrinos) of the T2K and NOz/A experiments, without relying on future second 
off-axis detectors (placed at a shorter distance or at the second oscillation maximum) and 
future upgraded proton luminosities (and/or detection technologies providing almost perfect 
detection efficiencies). 

The disadvantage is that the configuration at 14 km off-axis implies a 30% loss in statistics 
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with respect to the default configuration at 12 km off-axis. Consequently, the best strategy 
would be the following: one could tune the {E)/L of the two experiments to be the same by 
changing the T2K off axis-angle, since the design of the T2K experiment could in principle 
allow for such a modification, keeping the NOz/A off-axis distance fixed to 12 km which 
corresponds to (E) = 2.0 GeV. In order to reproduce the same {E)/L than the corresponding 
one for the NOi/A experiment, the energy at T2K should be 0.73 GeV, which could be easily 
achieved by placing the detector 2° off-axis, for which (E) — 0.75 GeV. In principle, since 
the configuration at 2° off-axis would also imply a gain in statistics for the T2K experiment 
with respect to the one at 2.5° off-axis, one would expect an improvement in the sensitivity. 
However, we noticed that the ratio for the two solutions Eq. (3) increases as the energy 
decreases. Since the configuration of NOz^A detector located at 14 km off-axis and T2K 
located at 2.5° off-axis provides a lower (E) than the energy for the configuration obtained 
by placing the NOi/A detector to its default location and moving slightly T2K to 2° off-axis, 
the sensitivity achieved in the former is higher than in the latter. 

The data fits we perform here are based on a analysis, where the function reads: 



where A'^"''* is the measured number of events, A'^*'*(sin^ 2^13, 5) is the expected number of 
events for a particular choice of the hierarchy and a refers to the statistical error correspond- 
ing to A^*'^. Our analysis includes the intrinsic Ue background. In our plots we show the 90% 
confidence level regions in the (sin^ 2^13, 5) plane resulting from the data fits, where the 90% 
CL is defined as |x^j„(wrong hierarchy) - x^^^(true hierarchy)] > 4.61 for 2 d.o.f [56]. It is 
very important to notice here that our 90% CL contours translate into the 97% CL contours 
for 1 d.o.f. ^ For all the sensitivity (exclusion) plots shown in this paper, the systematic 
error was not included but later we will make some comments on the impact of this error 
on our results. 

We summarize the results in Figs. 4, where we present the exclusion plots in the 

^ In order to test whether or not the 2 d.o.f approach is more appropriate, we use a Monte Carlo technique 
in which we simulate 10^ experiments for each point in the (sin^ 2^13, 6) parameter space and we compute 
the probability of distinguishing the hierarchy for each point. The result of the Monte Carlo analysis 
agrees with the conservative approach of 2 d.o.f but it does not reproduce the 1 d.o.f confidence levels. 
This Monte Carlo analysis has therefore guided us to use 2 d.o.f. 
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(sin^ 2^13, 5) plane for a measurement of the hierarchy at the 90% CL for the several pos- 
sible combinations, assuming that nature's solution is the normal hierarchy and Am|i = 
2.4 X 10-3 eV^ (left panel) and Am§i = 3.0 x 10"^ eV^ (right panel) (in the light of the 
recent MINOS results, we explore here the impact of a larger Amg-,^). Wc show as well the 
corresponding CHOOZ bound for sin^2^i3. A larger value of Amli implies more statistics 
and, consequently, a sensitivity improvement: see Fig. 4 (b), where for the sake of illustration 
only three representative configurations are shown. 

If both T2K and NOi^A run in their default configurations the combination of their future 
Phase I data (only neutrinos) will not contribute much to our knowledge of the neutrino 
hierarchy, sec the solid blue line in Figs. 4. If wc fix the T2K off-axis location to its default 
value of 2.5° but we change the location of the NOz/A detector to 14 km the improvement 
is really remarkable, see the dashed red line in Figs. 4: the sensitivity to the mass hierarchy 
has a milder dependence on the CP-phase S once that the (E) /L of the two experiments is 
chosen to be the same. The best sensitivity to the hierarchy extraction is clearly achieved 
when the NOi^A experiment is at 14 km off-axis and the T2K off-axis angle is the default 
one. If the T2K off-axis angle is slightly modified to 2°, see the dotted lines in Figs. 4 
it would be possible to reproduce the results from the combination of the data from T2K 
located at 2.5° off-axis and the NOi/A detector placed at 13 km off-axis. 

If the nature's choice for the neutrino mass ordering is the inverted hierarchy, the sensi- 
tivity curves depicted in Fig. 4 (a) will be shifted but in the opposite direction, see Fig. 5 
for the result from the combination of data from the Phase I of both experiments. Notice 
that the method presented here is independent on the octant in which 623 lies (if 6*23 is not 
maximal). There exists though a minor impact on the sensitivity curves since 623 > 7r/4 
(< 7r/4) implies more (less) statistics, see the first term in Eq. (2). 

We have explored the impact of an overall 4% systematic error. The sensitivity curves, 
while shifted slightly for the combination of NOi/A at 14 km off-axis and T2K at 2.5° off- 
axis (or NOi/A at 12 km and T2K at 2.0°), get much worst if both experiments run in their 
default configurations. One can understand this by making use of the bi-event neutrino- 
neutrino plots, see Figs. 2: a larger error bar has an enormous impact if the ellipses width 
is not negligible, sec Figs. 2 (a). 

The combination of data from an upgraded phase of the T2K and/or NOi/A experiments 
(by increasing the proton luminosities, the years of neutrino running and/or the mass of 
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FIG. 4: (a) 90% CL (2 d.o.f) hierarchy resolution for different possible combinations: the default 
one (T2K at an off-axis angle of 2.5° and NOuA far detector at 12 km off-axis, in solid blue), T2K 
at an off-axis angle of 2.5° and NOvA far detector at 13 km off-axis (long dash-dot red curve), 
T2K at an off-axis angle of 2.5° and NOvA far detector at 14 km off-axis (short dashed red curve), 
T2K at an off-axis angle of 2.5° and NOvA far detector at 16 km off-axis (three dots-three dashes 
blue curve) and T2K at an off-axis angle of 2° and NOvA far detector at 12 km off-axis (dotted 
blue curve). We have considered the statistics corresponding to the Phase I of both experiments. 
The vertical dashed line indicates the 95% CL CHOOZ bound for the value of Am'^^ for the panel, 
(b) The same as (a) but assuming that Am|-^ = 3.0 x 10^"^ eV"^ and only for the three most 
representative combinations: the default one (in solid blue), T2K at an off-axis angle of 2° and 
NOvA far detector at 12 km off-axis (dotted blue curve) and the optimal one, that is, T2K at an 
off-axis angle of 2.5° and NOvA far detector at 14 km off-axis (short dashed red curve). If one 
reinterprets these limits for 1 d.o.f then they correspond to the 95% CL, approximately. 

the far detectors) will obviously increase the statistics and will shift the sensitivity curves 
depicted in Fig. 4 (a), similarly to the effect of increasing Am-g^, see Figs. 6, where we have 
upgraded NOz/A and T2K experiments by increasing a factor of five their expected Phase I 
statistics (Phase II) . Fig. 7 depicts the results from an upgraded Phase II of both experiments 
in the inverted hierarchy nature's choice: if the neutrino mass hierarchy is inverted, the case 
for the Phase II of both experiments will be stronger, especially for Am|]^ = 2.4 x 10^^ eV^. 
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(a) Am^i = -2.4 x lO'^ eV^ 

FIG. 5: Same as Figs. 4 but assuming that the nature's choice for the neutrino mass spectrum is 
the inverted hierarchy. 
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FIG. 6: Same as Figs. 4 but increasing the statistics of T2K and NOuA by a factor of five. 

V. CONCLUSIONS 

The most promising way to extract the neutrino mass hierarchy is to make use of the 
matter effects in neutrino oscillations. For that purpose, the fastest way would be to exploit 
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CHOOZ 
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sin^ 2013 sin^ 2613 

(a) Am|i = -2.4 x lO'^ eV^ (b) Am^^ = -3.0 x lO'^ eV^ 

FIG. 7: Same as Figs. 6 hut assuming that the nature's choice for the neutrino mass spectrum is 
the inverted hierarchy. 

the neutrino data only from two near-term long baseline z/g appearance experiments per- 
formed at the same {E)/L, provided sin^2^i3 is within their sensitivity range or within the 
sensitivity range of the next-generation z/g disappearance reactor neutrino experiments. Such 
a possibility could be provided by the combination of the data from the Phase I of the T2K 
and NOi/A experiments. We conclude that the optimal configuration for these experiments 
would be 14 km off-axis for the NOi/A far detector and 2.5° off-axis for the T2K experiment. 
The combination of their expected results could provide a 90% confidence level (using 2 d.o.f ) 
resolution of the neutrino mass hierarchy if sin2 2^i3 > 0.11 (for Amj^ = 2.4 x lO'^ eV^) 
or if sin^26'i3 > 0.07 (for Aml^ = 3.0 x 10^^ eV^). A modest upgraded next Phase of both 
NOi/A and T2K experiments (by increasing a factor of five their expected Phase I statistics) 
could shift the 90% CL limits quoted above to sin^2^i3 > 0.03 (for Amj^ = 2.4 x lO'^ eV^) 
and to sin^ 26*13 > 0.025 (for Am|^ = 3.0 x 10"^ eV^). A slightly less sensitive combination 
is T2K at 2° off-axis angle and NOz/A at 12 km off-axis location. 
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